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A MANY-BODY RAGE THEOREM 


JONAS LAMPART AND MATHIEU LEWIN 


Abstract. We prove a generalized version of the RAGE theorem for 
JV-body quantum systems. The result states that only bound states 
of systems with 0 < n ^ N particles persist in the long time average. 
The limit is formulated by means of an appropriate weak topology for 
many-body systems, which was introduced by the second author in a 
previous work, and is based on reduced density matrices. This topology 
is connected to the weak-* topology of states on the algebras of canonical 
commutation or anti-commutation relations, and we give a formulation 
of our main result in this setting. 


1. Introduction and main result 

The RAGE theorem, due to Ruelle K 2 , Amrein- Georgescu [ 2 j andEnss [6] 
is a famous result relating the long time behavior of solutions to the Schro- 
dinger equation and the spectral properties of the corresponding self-adjoint 
Hamiltonian. In particular, it states that for any fixed x in the ambient 
Hilbert space Sj, 

1 f T 

lim — / \\KP c e~ ltH x \\ 2 dt = 0 , (1) 

T^oo T J 0 11 11 K J 

where P c is the spectral projector to the continuous spectral subspace of H, 
K is any compact operator and x(t) = e~ ltH x is the unique (weak) solution 
to Schrodinger’s equation 

j i x(t) = Hx(t), 

|x(0) = x. 

An equivalent way of formulating the same result is as follows: for every 
positive self-adjoint operator 7 in the trace class 6 * 1 (ij), consider the ergodic 
mean 

M(T) := 2 f e~ itH 'ye itH dt 
T Jo 

which is uniformly bounded in (S 1 ^). Then, the projection P C M(T)P C 
converges weakly-* to 0 as T A 00. Here, the operator 7 (t) = e ~ ltH ^ e ltH i s 
the unique (weak) solution to von Neumann’s formulation of Schrodinger’s 
equation 

J*7(*) = / 9 n 

\ 7 (0)=7- 

The previous formulation (PQ) corresponds to 7 = |x)(x|. Put differently, 
any weakly-* convergent subsequence of M(T ) has a limit M^ which is 
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supported on the space spanned by the eigenvectors of H. It can also be 
proved that commutes with H , that is, 


Mo o — CHj | ‘Pj ) (Pj | 

3 

where (pj is an orthonormal system of eigenvectors of H. 

The RAGE theorem is a very important result in quantum mechanics. For 
an infinite-dimensional Hamiltonian system such as Schrodinger’s equation, 
strong convergence to stationary states (that is, eigenvectors of H) cannot 
hold in general, due to the conservation laws. The RAGE theorem states 
that, on the contrary, weak convergence towards stationary states holds. 
Hence, in this sense, only bound states persist in the long time average. By 
virtue of its generality, self-adjointness of H being the only hypothesis, the 
RAGE theorem is a fundamental tool in the spectral theory of self-adjoint 
operators. For instance, it may serve as a first step towards a more precise 
scattering theory of quantum systems [ 16 ] and it is also often used in the 
study of Anderson localization mm- 

However, the information it provides for an interacting many-body system 
is often not very precise. Consider for instance three electrons in the field 
of a proton, described by the three-body Hamiltonian 


h 3 


( Ax i 


3= 1 



+ E 


1 

Xj -Xf\' 


( 3 ) 


It is known that H 3 has no bound state m flh] . The RAGE theorem 
therefore tells us that the ergodic mean M(T) tends weakly-* to zero for 
every initial condition 7 (a trace-class operator on T 2 (M 9 )). On the other 
hand, the corresponding Hamiltonian H2 for two electrons has finitely many 
bound states mmm and, of course, the hydrogen atom H\ for one electron 
has infinitely many. The physical picture is that some of the three particles 
escape, whereas the rest remain in a bound state of H-2 or H \, a phenomenon 
that is not captured at all by the RAGE theorem. The precise description of 
this process through scattering theory has been the object of several works 
m\IM- ^ was proved that, asymptotically, the time-evolved wavefunction 
can be approximated by a sum of products of the form < 8 > p where <p 

is a A-eigenfunction of H\ or H2 and 'tpt 0 . A tensor product of this form 
tends weakly to zero in L 2 (M 9 ), which is why the weak-* limit of M(T ) is 
always zero. 

In this paper, we would like to prove a new general version of the RAGE 
theorem that retains some information on the particles that do not escape, 
without addressing all the details of the scattering process. This is done 
by using another weak topology, for which the tensor product 7/y £3 p> con¬ 
verges to <p, and which was introduced by the second author in m- In this 
topology, the ergodic mean will converge to states with possibly less parti¬ 
cles, which are supported in the point spectrum of the Hamiltonians H n for 
n ^ N. This captures the principal physical ideas, even in situations where 
scattering is not known or not believed to hold, for instance for potentials 
with an arbitrarily slow decay at infinity. 
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Geometric convergence. In order to state our main result, let us quickly 
describe the notion of convergence used in [T 3 ], where all the details may 
be found. Let S) be any separable Hilbert space and denote by fj ;V the 
symmetric (or antisymmetric) IV-fold tensor product, respectively denoted 
by <8> s / a , of fj. Let Tfc be a sequence of IV-particle states, that is, r k ^ 
0 and Tr(rfc) = 1 . For instance, for a pure state F^ = |'Lfc)(’I / fc| for a 
normalized T/- € Sj N . The n-particle density matrix of F/,. is obtained by 
taking the partial trace with respect to N — n variables and multiplying by 
an appropriate normalization constant: 

tf ’ = (0) TW..MU)- 

The sequence T^ is said to converge geometrically to a state T if the reduced 
density matrices of all converge weakly-* to those of T. Except if con¬ 
vergence holds in trace-norm, the state T can never be an IV-particle state. 
It is necessary to work with states on the truncated Fock space 

F^ N (Sj) := C ©.£)©••• © S) N 

and this corresponds to the picture that some particles can be lost. For 
simplicity, all the states we consider in this paper are assumed to commute 
with the particle number, the theory for the general case is essentially the 
same. Such states can be written in block form as T = Go © • • • © Gn, and 
if we start with a sequence of IV-body states and investigate its geometric 
limits, these are the only states that can be obtained. 

We rephrase the previous discussion in the following: 

Definition 1 (Geometric convergence). A sequence of states {r / t}^ =1 on 
S) N converges geometrically to a state F = Go © ■ ■ ■ © Gn on F^ N (fj), if 



weakly-* in the trace-class © 1 (i 5 n ) for all n = 0 ,..., N. That is, Tr(KT^) —> 
Tr(AT(”)) for every compact operator K on Sj n . This notion of convergence 
is denoted as —*• g V and extended by linearity to sequences of states on 
the truncated Fock space F^ N (Sj). 

We emphasize that, by definition, the geometric limit T must always be 
a state. That is, it has to satisfy T ^ 0 and 

N 

Tr^v(r) = Go + ^2 ^"(^n) = 1- 

n= 1 

If convergence does not hold in then the final state has to live 

over spaces with less particles, but its trace is always equal to one. If all 
the particles are lost, then F is the vacuum state r = l©0©---©0. It is 
proved in m Lemma 3 ] that every sequence of states T has a geometrically 
convergent subsequence (the limit being a state). 
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Main result. We are now able to state our main result. We consider an 
abstract many-body Hamiltonian of the form 

n 

H n = ^ ^ hj + ^ ' wjg (4) 

j 1 

acting on the n-particle space S) n (for us symmetric or antisymmetric). Here 
h is a given self-adjoint operator acting on the one-particle space and hj 
acts on the jth factor. On the other hand, w is a self-adjoint operator on 
the two particle space f) 2 and Wji acts on the jth and tth factors. We make 
rather general assumptions on w in order to give a proper meaning to H n 
for all n. As will be clear from the rest of the paper, these can be weakened 


in specific examples or if one is only interested in a particular n = N. We 
assume that 

h is bounded from below (without loss of generality h > 1 ) ( 5 ) 

and that 

\w\ is infinitesimally form bounded with respect to h\ + h 2 (6) 
which means that 

s(hi + h 2 ) - C £ ^ w ^ e(/ii + h 2 ) + C £ (7) 


for all e > 0 . Under these assumptions, H n is bounded from below and may 
be realized as a self-adjoint operator for all n, by the method of Friedrichs. 

Theorem 1 (Many-body RAGE). In addition to ([ 5 ]) and ©, assume that 
for every compact operator K on Sj 

( K ! + K 2 )(h! + h 2 )~ 1/2 W (hi + h 2 )- 1/2 (8) 

is compact on S ) 2 , where K\ = K <g) and K 2 = 1^ <g> K. Let T ^ 0 with 
Tr(r) = 1 be a state on the (symmetric or antisymmetric) N-particle space 
Sj N . Then, the ergodic mean 

e -itH NTe itH N dt 

has geometrically convergent subsequences as T —> 00, and, for every such 
sequence, the limit is a convex combination of projections to eigenspaces of 
the n-body Hamiltonians H n for 0 ^ n ^ N. 

If we consider a sequence of times —>• 00 for which 

— / e-UHNYjiHN dt Mqo = Gq 0 ... 0 GjVj 

Tk Jo g 

the result states that 

G n — ^ a n,j\Vn,j)(Vn,j\ ( 9 ) 

j 

where {<p n ,j}j^i is an orthonormal system of eigenvectors of the n-body 
operator H n . By dehnition, geometric convergence means that the density 
matrices converge 
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weakly-* in the trace-class for all n = 0 ,...,N. Let us emphasize that, 

(n) 

because of the partial traces, the M > 0 are in general not supported on the 

(n) 

point spectrum of H n . The density matrix M > 0 should not be confused 
with the restriction G n of the state to the n-particle space. 

The condition © on w ensures that the interaction between a particle 
that stays and a particle that escapes vanishes in a weak sense. A condition 
of this type is clearly necessary to be sure that the remaining n particles are 
described by the Hamiltonian H n when the other N — n escape. The opera¬ 
tor B := (hi + h,2)- l / 2 w(hi + /12) -1 / 2 is always bounded by Assumption (|6|). 
However, in the applications, it is usually not compact. The picture is that 
w only decays in one direction, like the relative coordinate of the two parti¬ 
cles, and must be multiplied by a compact operator K in another variable, 
corresponding to the position of only one of the particles, as in (|S]) to make 
(K <g> 1 )B and (1 <g) K)B compact. 

Our theorem is stated for any initial datum T, possibly with an infinite 
energy. However, our proof does use the conservation of energy for smooth 
initial data. It is an interesting open problem to derive a similar result when 
h is not bounded from below. We will make more comments on this below. 

We also remark that the exact same theorem holds if the ergodic mean is 
replaced by 


1 

T 



e -iH N t Te iH N t dt 


where % is any nonnegative function such that f R y(f) dt = 1. 


Reformulation in terms of the CAR and CCR algebras. Theorem [Tj 
has a natural extension to the Fock space and the associated algebra of 
canonical (anti-) commutation relations (the latter is to be understood in the 
sense of Weyl operators, see [ 5 ] for a detailed introduction). This extension is 
based on the observation (see m Remark 6]) that for the rank-one operator 

K = |/l ®s/ a ■ ■ ■ ®s/a fn)(gi ®s/a ‘ ' ©s/a 9n\ , (10) 

geometric convergence of T*. G ©i(J r ^ JV ) just means that 

Tr = (gi ® s/a ■ ■ ■ ® s/a g n , r^/i <8> s / a ■ ■ ■ <8> s / a /n) 

= Trjr (a*(/i) • • • a*(f n )a(g 1 ) ■ ■ ■ a(g n )T k ) (11) 

converges. Here, a*(/),a(/) denote, respectively, bosonic or fermionic cre¬ 
ation and annihilation operators, depending on whether <8> s / a is the sym¬ 
metric or antisymmetric tensor product.We thus have: 


Corollary 2. Let A € (CCR(f)), CAR(ij)} be the C*-subalgebra of the 
bounded operators on J r (S)) satisfying canonical commutation relations, if 
is the symmetric Fock-space, or anti-commutation relations, if F{Sf) 
is the anti-symmetric Fock-space. Denote by Rv the projection of F{$)) to 
S) N and let p £ A! be a normal state on A, that is p(A) = Tr(T p A) for some 
Tp G ©1(7-”), satisfying [T p , n?v] = 0 for all A G N. Define the Hamiltonian 
on the Fock space F(S)) by 

H = 0 © h 0 Hjy. 

A/>2 
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Then the ergodic mean 

MA) :=^j\^ m Ae- im ) dt. 

has weakly-* convergent subsequences as T —>• oo and, for every such se¬ 
quence, the limit state is normal. For the density operator of this state, 
njvrooIlAr is a convex combination of projections to eigenspaces of for 
all N ^ 1 . 

Proof. The state px is normal for all T > 0 , with density matrix 

e~ im T p e mt dt = 0 i [ e~ iHNt U N T p U N e iHNt dt. 

Np 0 T J ° 

This sequence of normal states satisfies 

Tr^J © Y [ T e~ iHNt U N T p U N e iHNt dt 1 = £ T^ N {U N T p U N ) 

[a>AT 0 1 ) N^N 0 

which is small, uniformly in T, when Nq is large. By an ‘e /2 argument’, it 
is thus sufficient to prove the corollary for T p € For states on 

the truncated Fock space J~^ N , the weak-* convergence in A' is equivalent 
to geometric convergence m Remark 6]. In the fermionic case this is an 
immediate consequence of the formula dill) and the commutation relations. 
In the bosonic case one uses additionally that on S) N the Weyl-operators 
may be expressed as convergent power series in creation and annihilation 
operators (see [ 5 j chapter 9 ]). Therefore, the existence of weak-* convergent 
subsequences as well as the properties of their limits follow from Theorem [lj 
by linearity. □ 

Note that A = CCR(ij) is not separable, hence bounded sequences do not 
necessarily have weakly-* convergent subsequences. However, we are able to 
prove that our particular sequence px has convergent subsequences, using 
the fact that it is obtained from a fixed normal state by a particle-number 
conserving evolution. 

Application to Schrodinger operators. A typical example to which our 
result applies is that of a non-relativistic system of N fermions or bosons, for 
which the interaction w is a function of the relative position of the particles, 
that tends to zero at infinity in a weak sense. 

Corollary 3 (Non-relativistic Schrodinger operators). Let Sj = L 2 (R d ), 
h = — A + R(x) + e and w be the multiplication operator by an even function 
w(x i — X2). We assume that V = fi-\- f2 cmd that w = fj, + f± where 
fi G L Pi (R d ) for some max(l,d/ 2 ) < pi < 00, or fi e L°°(R d ) and fi 0 
at infinity. Then Theorem\T\ holds for the many-body Hamiltonian 

N 

h n = y j ~A Xj +v(xj)+ y w ( x j - x i ) ■ 

j =1 l^j<£^N 
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Proof. Under the assumptions of the corollary, (1 — A) _1 / 2 U(x)(l — A) -1 / 2 
and (1 — A) _1 / 2 ro(x)(l — A) -1 / 2 are compact on L 2 (M rf ) [31 Chap 8], hence 
the hypothesis ©, ([6|) are satisfied and h > 1 for an appropriate choice of 

e. 

The verification of ([8]) is more involved, the intuition however is rather 
simple. It relies on the fact that (hi + h2)~ l ^ 2 w(hi + is compact in 

the relative coordinate x\ — X2 parametrizing the subspace V = {(x\,X2) € 
M? d : x\ + X2 = 0 }, i.e. its action on the space L 2 (V) is compact. This, 
together with the fact that K 0 1 is compact in the direction W l x { 0 } 
parametrized by x\ and (M rf x { 0 }) 0 V = M 2d , implies that the product is 
compact. More precisely, the property we use is the content of the following 
lemma, which we prove in Appendix lAl 

Lemma 4 . There exist compact operators Kj on L 2 (V) and bounded oper¬ 
ators Bj on L 2 (V 2 ~) such that 

OO 

(h\ + h 2 )~^ 2 w(hi + h 2 )' 1 ^ 2 = ^2 Kj 0y Bj (12) 

3 = 1 

and the sum converges in the operator-norm. 

The notation 0y emphasizes that the tensor product of operators is in¬ 
duced by L 2 (R 2d ) = L 2 (V) 0 LPiy 2 -). The property (fT 2 h is known as V- 
semicompactness and implies hypothesis © by [U Proposition 9 . 2 . 2 ], 

The argument is as follows: It is clearly sufficient to prove that (K 0 
1 )(A 'j 0y 1 ) is compact for every j. Then, since Hilbert-Schmidt operators 
are dense in the compact operators, it suffices to show that C := (A 0 
1 )(B 0y 1 ) is compact (and actually Hilbert-Schmidt) for operators A € 
©2(-h 2 (®l d )), B € & 2 (B 2 (V)). To prove this, let a(x,x'), b(v,v') be the 
integral kernels of A and B , respectively. The operator C acts of if € 
L 2 ( R^xK^as 

(Cij})(xi,x 2 )= / a(xi,x')h(x' — X2,x' — y')'if(x , ,y') dx'dy'. 

J E 2d 

So C is an integral operator, whose kernel is easily seen to be in L 2 (R 4rf ). □ 


If w and V decay fast enough at infinity, a precise theory of scattering 
is available pa Hugos]. The existence of the wave operators and their 
completeness implies that there exists vectors 'i/Vtj such that 


lim 

t —^-)-oo 


D -itH N 


N 

71=0 j 


g I 


itA , 


-n) 


® ^Pn,j 


= 0 


L 2 ( R dN ) 


( 13 ) 

where H n (p n j = A n ,j<Pn,j- For the ergodic mean of the time-evolved projec¬ 
tion |'h)( 4 '| we obtain 


1 

T 



N 


e ~itH N 


|T)('I'| e itHN 



a n,j I Pn,j ) (Pn,j \ 
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geometrically, where a n j = ||^ n j || 2 gives the coefficients in Q. Additionally, 
any combination of the eigenvectors tpj jn of the n-particle Hamiltonians H n 
can occur in the limit, so any Go,..., Gat of the form ([9]) may be obtained 
in the geometric limit, by choosing an appropriate initial condition. We do 
not know if the same property holds under the more general assumptions of 
Theorem |T] 

Results similar to Corollary [3] hold for Schrodinger operators with mag¬ 
netic fields, for pseudo-relativistic operators, etc. We do not state them here 
for shortness. 

The rest of the paper will be devoted to the proof of Theorem [TJ 


2. Proof of Theorem □ 

Let T be a state on the IV-particle space S) N and T k —^ oo be a sequence 
of times such that the ergodic means 

1 r^k 

M(T k ) := J e~ itHN Te itHN dt 

with initial condition T converge geometrically to a limit M a0 = ©^ = 0 G n . 
Such a sequence exists for every T by m Lemma 3]. 

The proof of Theorem [T| proceeds by showing that M 0 0 is left invariant by 
the Hamiltonian HI := ®n=o^n (with Hq = 0 and H\ = h) on the truncated 
Fock space = C ®^ =] ij n . That is, we have M 00 = e~ zm M 00 e im and 
thus e~ ltHn G n e ltHn = G n , for every t € M. Consequently, the eigenspaces of 
G n are H n -invariant. As G n is trace-class, the eigenspace corresponding to 
a non-zero eigenvalue of G n has finite dimension, and is thus a direct sum 
of eigenspaces of H n . 


Step 1. The first step to proving invariance of M 0 0 is to note that, for every 
s € M, the sequence e~ lsHN M(T k )e lsHN also converges to M 0 0 geometrically. 
This holds because geometric convergence is controlled by the trace-norm 
and 


1 / k 

e -isH NM ^ Tk ysH N = _ / e -i(t+s)H NTe i(t+s)H N dt 

T k Jo 

e -itH NTe itH N dt 

1 / /'T i + S r s 

= M(T k ) + — I / e~ itHN Te itHN dt- e~ itHN Te itHN dt 

T k \Jr k Jo 

^ - 

^2 s/T k in Si 



g 


(14) 


Having established this, we would like to prove that also 

e - is H NM (T k ) e is H N e - im Mooe ism . (15) 

g 


As the left hand side depends only on H jy, which equals the restriction of 
HI to Sj N , and the right may depend on HI on all the sectors with n ^ N 
particles, this will certainly not be true for arbitrary Hamiltonians on Fock 
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space. In fact, the proof of (fT5l) will depend crucially both on the properties 
of H and the sequence M(T]f). 

Step 2. Using the particular form of the sequence M (7).) we can reduce 
the proof of (fl5l) to sequences of bounded energy. That is, let be such 
that Hn + ejv > 1 and assume the initial condition M(0) = T satisfies 

TV ^H n + e N f/ 2 T{H N + e N ) 1 / 2 ) < E 

for some constant E > 0. Then, since e~ ltHN commutes with Hn and 
preserves the trace-norm, we have 

TV ({H n + e N ) l ' 2 M{T k ){H N + e N ) l > 2 ) < E , (16) 

for every k € N. Since any initial condition M(0) = T can be approximated 
to arbitrary precision by states of finite energy, it is sufficient to prove (1151) 
for sequences satisfying (1161) . by virtue of the following lemma. 


Lemma 5. There exists a constant c, depending only on N, such that 
for every sequence { 7 k}kLi in © 1 (iy Y ) that converges geometrically to 7 € 
&i(E^ n ) and satisfies Ibfclleu^) < £ we have IMIe^^JV) < ce. 


Proof. Since the trace-class is the dual of the compact operators, we have 


7 


(n) 


= sup lim 

®i(h”) incompact, |if|=1 k—Kx> 


TV? 


fix 


(n) 

k 



sup lim |TV fl jv (K ® l A r_ n 7 fc )| ^ 

|r:||=i fe->oo 



e. 


Now 7 = ®^_ 0 G n is completely determined by its reduced density matrices 
(cf. [T3J Lemma 1 ]), explicitly 


G n = 7 (n) +E(- 1 ) j 

3 = 1 


n + j 
n 


TV 


n+l,...,n+j 7 


■(n+j) 


which proves the claim. 


□ 


It will often be useful to state the energy bound using a Hamiltonian 
without interaction. Let 

n 

3 = 1 

and 

■= {t € ©i(J5") : \\Vh*tVh?\\ Gi < . 

This space has a natural norm given by I 7 I := | 7 || s + || H^\\q 1 . 

Using H n instead of H® gives an equivalent norm, because of the inequali¬ 
ties ©• We also define 

©i.ho^) == {7 € < 5 i(E <N ) : || < 00 } , 

with H 0 = Hq = 0. It will be important that geometric conver¬ 

gence preserves such energy estimates. The following lemma proves that the 
unit ball of Snof^) is closed under geometric convergence. This is the 
only step in the proof of Theorem [T] for which positivity of H® is essential. 
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Lemma 6. Let {p'k}kL\ i n ©1 u°(^ N ) be a bounded sequence that converges 
geometrically to 7 . Then for every n ^ N 



in Si(ij n ). 


Proof. By linearity, it is sufficient to prove the claim for 7 ^ £ &ih° 
For any n ^ N we have 

0 ^ (H° n ® v - n)^)- 1 = 1 - (v ® K- JiKr 1 < 1, 

because H® ^ 1. Thus, for an arbitrary g € &i H o(S) n ), 



So \J #° 7 fc yjH® is uniformly bounded in ©i(f) n ) and there exists a sub¬ 
sequence such that 



Now for every compact operator K on Sj n we have 


lAyW) = lim Tr (H*) 1 / 2 ) 

K —>00 

= Tr (A7 (n) ), 


whence 7 ^ H^ n ' ) y/H®. 


□ 


Step 3. The reduction of the problem to sequences of bounded energy in 
step 2 will now allow us to study the sequence e~ lsHN M{T}f)e lsHN via the dif¬ 
ferential equation it satisfies. If 7 (f) is a solution to the von Neumann equa¬ 
tion ([2D, its reduced density matrices (formally) satisfy the finite BBGKY 
hierarchy 




[#n,7 (n) W] + (n + l)J]Tr n+ i (ju^n+i,7 (n+1) (*)]) • 

j =1 


Note, however, that the equation has no clear meaning (not even in a weak 
sense) if w is not bounded, due to the partial trace. If w is a bounded 
operator D((Ho) a ) —» ij 2 for some a > 0 , this problem can be handled by 
considering only initial conditions satisfying Ti-((H^ r ) a 'y(H ( ^) a ) < E, but to 
deal with potentials that are only i^-form-bounded we will have to define 
a modified equation, that is equivalent to the original one for bounded w. 


Proposition 7 (Well-posedness of the truncated BBGKY hierarchy). 
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( 1 ) For every 7 € 6 1H o(J r ^ Ar ) the family of reduced density matri¬ 
ces 7 i s the unique solution {7 ( n \t) : n = 
0 ,... ,N} to the the system of equations 


7 ( n )(t) =e -itH„^(n) e itH„ 


n pt 

- *C n + X) E J 0 e- i{t ~ s)Hn Tr „ +1 (c jn (w^ n+l \s))) e^^ds, 


£jn (w,7 ( ” +1) ) = f^n+l 2w J,n+l^n+l 2 .^+l7 (n+1) (s)^i+l 


( 18 ) 


such that 


7 W(^r(R, 6 1 , H o(^)). 

( 2 ) For erer?/ bounded sequence 'y k € 0 li uo(J rsSAr ) that converges geomet¬ 
rically to 700, f/ie corresponding solutions 7fc(i) of the von Neumann 
equation © with initial condition 7^ converge geometrically to the 
solution with initial condition 700/ 


Vf <E M : 


e~ im ^ k e im 


_a e _ * ffl 7ooe* ffl . 
g 


Proof. Note that 7( Ar +V = 0 and this is a triangular system, which can be 
solved starting with 7 ( N \t) = e~ ltHN ^ N ) e ltHN . This immediately gives us 
uniqueness, for assume we have two solutions with 7^'^ (0) = 7W = 0^(0), 
then 7 W(t) = 0( N )(t). If € L°°(M, 6 1 iH o) this then implies 

7 (^- 1 ) (t) = 0 (^- 1 ) (t). 

We now check that 7(£)( n ) = j s indeed a solution. By 

linearity, we may restrict to initial conditions in 7 € ©1 h°(^> N )- Since both 
(H n + ejv ) 1/2 (^)- 1/2 and {H%y/ 2 {H N + e N )~ x ' 2 are bounded, due to 
hypothesis (J6]) , we have 7(t) = 7 (t)( N ) € L°°{R,& 1 H o(S) N )). Equation (fTTli 
then implies that 7 (t)^ € L°°(R, ^o(^ n )). 

Now assume for a moment that w is bounded and let 

W n := ^2 W F 
1 


be the interaction of the first n particles with the remaining N — n. Then 
we have 

Hn = H n <8> l^iv-n + 1 fin <8> Hn-ti + W n , 
and we can write 7 (t)^ using Duhamel’s formula, 

7(t) (N) = g itHN—n g itH n git.HN — n 


-*/ 


e ~i(t—s)H N - n e ~i{t-s)H n 


W n , 7(s) (JV) l e dJs)H N -n e i{t-s)H nds 
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(with the 1 ( 8 > omitted for shortness). This implies that 

TC*)(”) ~ e~ itHnr y ( ' n ' > e itHn 


= — % 


N 


n 


e -i(t-s)H n rjy 


n+l,...,N 


rt n 


= —l 


i{n + 1) / E e -i(t-s)H n Tr 


n +1 


,0 i=l 


W n ,7(s) (A °J) e i{t ~ s)Hn ds 

V+1 , 7 (s) (ri+1) ])e^-^d S . 


As /qE7(s)("-+') hE, is trace-class and in is bounded, we have 


TWi 


Wj,n+ 1,7 



= Tm+i 


,- 1/2 

' l n +1 w j,n+l 



which gives (fTKl) . In order to account for unbounded interactions, let w s := 
(SH® + l)~ 1 ^ 2 w(6H2 +1) -1 / 2 for 5 > 0, which is bounded by hypothesis ([6]). 
We clearly have 

lim (H^-^w^Hg)- 1 / 2 = (H^-^wiH^)-^ 2 

< 5->0 


in the strong topology of operators on Sy. Hence, the constant C £ of equa¬ 
tion (ED may be chosen in such a way that the inequalities hold for all 
0 < 5 ^ 1 . Then, for an appropriate choice of ^ 0 , 

(<+ E ™% + e N y 1/2 ■■Z N ^D((H° n ) 1 / 2 ) 


is bounded, uniformly in 5 . This implies that, as operators on D^(H ^) 1 / 2 ), 

hm(<+ El w j£ + = (Hn + &n) 

1 


strongly. This in turn implies that the unitary groups Us generated by these 
operators also converge in the strong operator topology on D^H ^) 1 / 2 ). So, 
for 7 € 6] , we have 

lim Us (i) 7 = (H%) 1 / 2 e~ itHN 1 

< 5—>-0 

in trace-norm. We thus have, using (fT 71 ) , that 

\im(Us{t)^Us{t)*) {n) =7 (n) (*) 

< 5->0 

in &i As w s is bounded, the left hand side solves (fl8l) with this 
interaction. To take the limit on the right hand side of ( 1181 ) observe that 

lim Tr n+1 C jn ((u/ 5 , 7 (ri+1) )) = Tr n+1 Cj n (E, 7 (n+1) )) 


in the space of operators for which (H®) 1 / 2 A{H^) x / 2 is trace-class. As 
£7 ^0 is obviously contained in this space, this shows that 7 (t)^ is indeed 
a solution. 

We now prove the continuous dependence on the initial condition, item ( 2 ). 
Since we have already proved uniqueness, it is sufficient to show that the 
geometric limit of 7fc(f) = e~ itH 7 fc e * <H is a solution of ( 1181 ) . On the space 
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Sj N this is trivial, since for any compact operator K , K(t) := e itHN Ke iIHn 
is also compact and thus 


lira Tr (Ke~ itHN ^ k e itHN 

k—yoo 


) = Tr {Ke~ itHN ') 00 e itHN 


) 


for every t € R. Now let m < N and assume that —*■* 7 oo(t) (n) 

in ©i(f) n ) for every t € M and n = m + 1 ,... N. Then, we have for every 
compact operator K on S) m 

(tv( ir 7fe (0 (m) ) -Tisim ^e~ itHrn K^/k (0) (m ^ e itHm )) 

ft 771 

= / XI Tr b m+1 (( K ( f ~ S ) ®s/a 1 )£ J - m (w,7 (m+1) )) ds. 

• y ° 1=1 

It is clearly sufficient to prove convergence for K in a dense set of compact 
operators, so we may assume that K(t—s ) = K(t—s){H m )~ 1 ^ 2 with compact 
K. For such an operator K , the integrand is uniformly bounded, and we will 
show that it converges pointwise. Using Lemma 0 the induction hypothesis 
gives us 

( hL .) 1 / 2 tL +i| (*)« +1 ) 1/2 =. « + 1 ) 1 / 2 t £“ +1| (»)« +1 ) i/2 . 

It is thus enough to prove that 

(kj{t - s) ® s/a l){h m+1 hjy 1/2 w jt m +1 {H^ +1 )- 1/2 , 

where Kj(t — s) = K(t — s)(H^) _1 / 2 /iV 2 , is compact. Hypothesis ((HJ) guar¬ 
antees that 


[L <S> s /a l)(^m+l + hj ) ^ ^ 

is compact for any compact operator L. Hence, it is sufficient to show that 

OO 

(Kj(t - s) <8) l)(h m+ ihj)~ 1/2 (hj + h m+ i) 1/2 = ^ L; <g> B t , 

i=i 

as a norm-convergent sum, with bounded operators B[ on f) and compact 
operators L[ on Sj m . Since hj and h m+ \ commute, we have 

((fc m +1 + M/fti'wi) -172 = (i + fc- 1 )" 1 ' 2 h - LT§1 

and because 0 < (1 — /i“ 1 )/(l + hj 1 ) < 1 we can write this using the 
convergent power series of (1 — x)^ 1 / 2 on |x| < 1, which gives the desired 
form. 

We have thus shown that 7converges weakly-* to the right hand 

side of equation (fT8l) with 7which proves the 
claim by the uniqueness of solutions to (fT8|) . □ 



This completes the proof of Theorem [TJ To summarize: We established 
that the Theorem is implied by invariance of the limit state under 
EL As a first step m, we then showed that this holds if the limit of 
the time-evolved sequence is the evolved limit state (equation (1X511 ). In 
the second step, we reduced the problem to considering initial conditions 
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with bounded energy, and thus to proving (fl5l) for sequences of uniformly 
bounded energy, i.e. satisfying (fT6l) . The final step consisted in studying 
the BBGKY-hierarchy satisfied by the reduced density matrices under the 
condition of finite energy, in its weak form (1181) . We proved uniqueness 
of solutions to this equation in Proposition [71 Additionally, we proved that 
that geometric convergence for the initial conditions, together with a uniform 
energy bound, implies geometric convergence of the corresponding solutions 
at any finite time. For the sequence M(T k ) of ergodic means this gives 

Or^-i © e~ iHNt M(T k )e iHNt = e~ m (0 © M(T k ))e mt , 

g 

that is (1151) . As we have shown, this implies invariance of the limit state 

e-^Mooe*® = g-lim M{T k )e iHNt ® , 

and thus Theorem [I] 

Remark 1. Instead of studying the evolution equation for the density ma¬ 
trices as in step 3, one could also study the dual evolution on the space of 
operators that we test against. The condition that needs to be verified is that 
this space is invariant under conjugation by 

Geometric convergence of 7 k is defined via convergence of Tr((iF <S> s / a 
l) 7 fc); so the natural space of operators is K <S> s / a 1, where K is compact on 
Sj n for some n ^ N. However, this space will in general not be invariant 
under e~ im (e.g. if w is not bounded, the time-derivative is never tangent 
to this space). This amounts to the fact that the map 7 i->- i s no t 

continuous with respect to the geometric topology on ©1 (P^ N ). However, 
Proposition [7] proves that this map is continuous for the induced topology on 
bounded sets of &i^o(P^ N ). These sets are closed by Lemma® The dual of 
this space is spanned by the operators A® s / a 1 such that (iL°) _ 1 / 2 A(iL °) _1 / 2 
is compact, and it can be shown to be invariant by a reasoning similar to 
that applied in part two of Proposition [?| 

This is sufficient to obtain Theorem [ 7 ] for arbitrary initial conditions, 
since the sequence M(T k ) can be uniformly approximated in that space. 

Appendix A. Proof of Lemma ® 

Recall the notation V := {(xi, X 2 ) € M. d X M d : x\ + x '2 = 0} and parame¬ 
trize this subspace by the coordinate v = x\ — X 2 - We will begin by showing 
that 

OO 

(1 - A X1 - \ X 2 )- l ' 2 w{v)(l - A ,)- 1 / 2 = Y, K 3 B 3 » (19) 

3 =1 

with compact I\j and bounded Bj. 

Let £ denote the conjugate Fourier-variable to v and ( that to r 1 := 
X 2 + x\. Let rj(£X) = (1 + 2 £ 2 + 2C 2 ) -1 ^ 2 , so (1 — A ) -1 / 2 is just the 
Fourier-multiplier by rj. Since 7 tends to zero at infinity, there exist functions 

0 with compact support, such that Y!j=i fjiOdjiO r l(f- C) in 
L°°. We can additionally arrange to have Ylj =1 fj(OSj(C) ^ v (€5 0 f° r every 
k. Now define Bj := Pf ] gj(f)J- v ± and Kj := P)) 1 fj(f)J r v w(v)( 1 —A^)~ 1//2 , 
where T x is the Fourier transform in the variable x. These operators clearly 
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have the desired properties, and we now need to check convergence of the 
sum CD, which is not immediately obvious if w is not bounded. First, note 
that Kq := J 7 ^ 1 ( 1 + F v w{y){ 1 — /S. v )~ l t 2 is a compact operator on 

L 2 {V ) and f n y/l+£ 2 is bounded. In order to exploit the fact that ATq y 1 
is compact in the first factor, we identify the operators on L 2 {V) <S> L 2 (V ri -) 
with operators from L 2 {V) to C(L 2 (V' L ),L 2 ( R 2d )) via A(ip)(p := Aifj <8) ip. 
After this identification, Kq ®y 1 dehnes a compact operator, since for any 
weakly convergent sequence ipk —*■ if) in L 2 (V ) and tp € iAiV 2 -) 

\(K® V l)(i/> k )<p - (K ® l)(i/;)<p\\ ^ WK't/jk - Ktp\\ L2{v) |^| L 2 (y ± ) • 

Furthermore, the operators 

m 

Am ■.= Y J ^~ l fAoV^+e9j{cv 

3 =1 


converge to T £)-^/l + £ 2 F in the strong operator topology of 
£(L 2 (V), £(L 2 (V 1 -), L 2 ( R 2d ))), since for every ip € L 2 (V) 


lim 

m—>-oo 


An(V0 - ^ C) V(^ 


= lim sup 

m—> oo . 


-v(€,0) v / i+^ 2 -0(O 


£(L 2 (R J -),I/ 2 (R 2li )) 

2 

^ = 0, 


1=1 


as the integrand is bounded by |2r/y'T+£ 2 V’| 2 € L 1 and converges to zero, 
pointwise in £ and uniformly in £. Thus, since Kj®yBj = A m (KQ®y 1), 
the sum in (11911 converges in norm. To deduce the statement of lemma [4] 
from (fT9l) we use use the formula 

(hi + h 2 )~ 1/2 - (2e - A X1 - A X2 ) -1 / 2 


( 20 ) 


J z 1/2 {hi + h 2 - z) 1 (14(aq) + V(x 2 )) (2e — A X1 — A X2 — z) 1 dz, 


where a is the boundary of a sector in the right half plane Re(z) > 0 with 
ad {Im( 2 ;) = 0} = {e} and the integral converges in the operator norm. The 
argument we used to prove (fl9l) also implies that (1 — A Xl — z)~ l / 2 V (aq)(2e— 
A X] — A X2 — z)^ 1 / 2 is x 1 -semicompact, and thus (cf. Lemma ED 

(1 - A X1 - zY l/2 V(xi){2e - A X1 - A X2 - 2 ) _1 w(u)(l - A„)~ 1/2 


is a compact operator for every z € a. Now for / € L p (M d ) with p > 
max(l, d/2) the Kato-Seiler-Simon inequality gives 


(l-A x -z) 1/2 /(x)(l - A x - z) 1/2 


^C( Pl d)(l+\z\ 2 ) W4p - 1/2) \\f\\ LP . 

6 P 


Using the decomposition V = f\ + f 2 we see that the integral in (12U1) still 
converges in norm if multiplied by iu(u)(l — A„)~ 1//2 from the right. Con¬ 
sequently, it defines a compact operator, so the difference between equa¬ 
tions (fT2ll and (fT9l) is compact and the proof is complete. 
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